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Abstract 

This paper presents a general framework for analytically obtaining the closed-form probability density function 
(PDF) of the Euclidean distance between any arbitrary reference point and its n-th neighbour node, when N 
nodes are uniformly and independently distributed inside a regular L-sided polygon. First, we exploit the rotational 
symmetry of the regular polygons and quantify the effect of polygon sides and vertices on the distance distributions. 
Then we propose an algorithm to determine the distance distributions given any arbitrary location of the reference 
point inside the polygon. For the special case when the arbitrary reference point is located at the center of 
the polygon, our framework reproduces the existing result in the literature. The proposed framework has many 
implications for the modeling and analysis of wireless networks. 

Index Terms 

Wireless networks, random distances, distance distributions, regular polygons. 

I. Introduction 

Recently, the distance distributions in wireless networks have received a lot of attention in the litera- 
ture [|T|-[[3|. The distance distributions can be applied to study important wireless network characteristics 
such as interference, outage probability, connectivity, routing and energy consumption [[T|, Q-Q- 
distance distributions in wireless networks are dependent on the location of the nodes, which are seen 
as realizations of some spatial point process. When the node locations follows an infinite homogeneous 
Poisson point process, the probability density function (PDF) of the Euclidean distance between a point 
and its n-th neighbour node follows a generalised Gamma distribution [|7|. However, as recently identified 
in p]], pj, [[8|, [[9|, this model does not accurately reflect the distance distributions in many practical 
wireless networks where a fixed and finite number of nodes are uniformly and independently distributed 
over a finite area such as square or hexagon or disk region. Note that these finite regions of interest can 
be conveniently modeled as special cases of a regular L-sided polygon (referred to as L-gon for brevity). 
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e.g., L = 3, 4, 6, oc correspond to equilateral triangle, square, hexagon and disk respectively. In this 
context, the two important distance distributions are: (i) the PDF of the Euclidean distance between two 
nodes uniformly and independently distributed inside a L-gon and (ii) the PDF of the Euclidean distance 
between any arbitrary reference point and its n-th neighbour node, when nodes are uniformly and 
independently distributed inside a L-gon. For the first case, the PDF of the distance between two nodes 
uniformly and independently distributed inside an equilateral triangle pO| |, square [ [TT| , [[12), hexagon [[3|, 



[ [T3] | and disk JTT] ] are well known in the literature. These results are, in fact, special cases of the general 



result obtained recently in [[14). For the second case, the PDF of the Euclidean distance to the n-th 
neighbour node is obtained in [fp] for the special case when the reference point is located at the center of 
the L-gon. 

In this correspondence, we present a general framework for analytically obtaining the closed-form PDF 
of the Euclidean distance between any arbitrary reference point inside the L-gon and its n-th neighbour 
node. The proposed framework is based on characterising the overlap area between the L-gon and a disk 
centered at the arbitrary reference point located inside the L-gon. There are two key insights which lead 
to our results: the use of the rotation operator that simplifies the characterisation of distances and overlap 
areas, and the systematic analysis of the effect of L-gon sides and vertices on the overlap area. Based 
on our proposed framework, we formulate an algorithm to determine the distance distributions given any 
location of the reference point inside the polygon. We provide examples to demonstrate the generality of 
our proposed framework. We also show that the result in [[T| can be obtained as a special case in our 
framework. 

The following notations will be used in the correspondence: Boldface upper and lower case letters 
denote matrices and vectors, respectively, [x, yY denotes the coordinates of a point where [-Y denotes 
transpose of a vector. [-J denotes integer floor function, abs(-) denotes the absolute value or modulus. 
min(-) denotes the minimum value and max(-) denotes the maximum value, aj denotes the j-th element 
of vector a. det(-) denotes the determinant, ^4=1^1 denotes the area of the two dimensional region A, 

II. System Model 

Consider nodes which are uniformly and independently distributed inside a regular L- sided polygon 
^ G M^, where denotes the two dimensional Euclidean domain. Let u = [x, yY ^ A denote an 
arbitrary reference point located inside the L-gon. 



A. Polygon geometry 

Without loss of generality, we assume that the L-gon is inscribed in a circle of radius R and is centered 
at the origin [0, 0]^. Then, its inradius is Ri = Rcos{ti/1) and its area A is given by 

A=\A\ = -LR^^m(^\ . (1) 



2 V ^ , 

Let Si and denotes the sides and vertices of the polygon, for ^ = 1, 2, . . . , L, which are numbered 
in anti-clockwise direction, as shown in Fig.[T] We assume that the first vertex Vi of the polygon is at 
[R, 0]^, i.e., at the intersection of the circle inscribing the polygon and the x-axis. The interior angle of 
the polygon 6 and the central angle between two adjacent vertices i9 are given by 

e=-J^ and tf=^. (2) 

B. Rotation operator 

For compact representation, we define the rotation operator which rotates an arbitrary point u = 
[x, yY anti-clockwise around the origin by an angle The rotated point Wu can thus be expressed as 

{^^u) = Tu, (3) 

where T is the corresponding rotation matrix of and is given by 

cos(M) -sin(M) \ 

. (4) 

sin(M) cos(M) I 

Also define 9^"^ as the inverse rotation operator with rotation matrix T"^, which rotates an arbitrary point 
u = [x, nY clockwise around the origin by an angle I'd. We note that the rotation under the operator 
is an isometric operation which preserves the distances between any two points in the two dimensional 
plane, i.e., det(T) = 1 and TT"^ = I, where I denotes the identity matrix. 

C. Distances from the polygon sides and vertices 

In this subsection, we find the distances from an arbitrary reference point u located inside the polygon 
to all its sides and vertices. First, we examine the distance to any vertex. Using the geometry, the distance 
d{u] Vi) between the point u and the vertex Vi is given by 

d{u', V^) = ^/{x-Rf + y\ (5) 



In order to find the distances to the remaining vertices, we use d{u] Vi) and exploit the rotational 
symmetry of the L-gon. By appropriately rotating the point u and then finding its distance from vertex 
Vi, we can in fact, find the distance to other vertices. Thus, using ([5]) and the rotation operator defined 
in ([3]), we can express d{u] Vi) as 

d{u] Vi) = d{dl-^^-^^u] Vi). (6) 

Next, we examine the distance to any side. We define this distance to any side as the shortest distance 
between the point u and the line segment formed by the side of the polygon. If the projection of the 
point u onto a side lies inside the side, then the shortest distance is given by the perpendicular distance 
between the point and the side. If the projection of the point u onto a side does not lie on the side, 
then the shortest distance is the minimum of the distance to the side endpoints. Using the geometry, the 
perpendicular distance from an arbitrary point u to the line segment formed by side Si is given by 

abs (y + tan (^) x — i?tan f^)) 
Sr) = , (7) 



^l + tan^d) 

and the distance to the side Si endpoints is simply given by d{u] Vi) and d{u; V2). Thus, the shortest 
distance d{u] Si) to the side Si can be expressed as 

{mm{d{u] Vi), d{u] V2)), max((i(i6i; l/i), d{ui]V2)) > t] 
(8) 
p{u] Si), otherwise; 

where t = 2i?sin(^) is the side length of the ^-gon. Then, using ([3]), we can express p{u] S^) and 
d{u] Si) as 

p{u; Si) = p{^-^'-'^u; Si), d{u; Si) = d{^-^'-'^u; Si), (9) 

Definition 1: Define the distance vector d as an indexed vector of the distances from the arbitrary point 
1^ G ^ to all the sides and vertices of the L-gon, defined in ([5]), ([6]), ([8]), and (|9]) as 

d = [d{u; Si), d{u; Sl)^ d{u; Vi), d{u; Vl)]. (10) 

III. Problem Formulation 

Consider a disk V{u;r) of radius r centered at the arbitrary reference point u. First, we define the 
probability that a random node, which is uniformly distributed inside the L-gon A, lies at a distance of 
less than or equal to r from the point u. 



Definition 2: Define the cumulative density function (CDF) F{u]r), which is the probabihty that a 
random node falls inside a disk V{u]r) of radius r centered at the arbitrary reference point u, as 

where 0{u] r) = \T){u] r) n A\ is the overlap area between the disk V{u] r) and L-gon A. 

Then the PDF fn{u]r) of the distance from an arbitrary point u to the n-th neighbor node is given 

by 

_ {l-F{u-r))''-{F{u-r)r-' dF{u-r) 
J^^""'"^^ - B{N-n + l,k) dr ^^^^ 

where is the number nodes are which are uniformly and independently distributed inside the L-gon, 
B{a,h) = r(a)r(6)/r(a + h) is the beta function, r(-) denotes the gamma function, ^(•) denotes the 
partial derivative respect to the variable r and F{u] r) is given in ( [TT] ). 

The challenge in evaluating ( pTj ) and ([12]) is finding the overlap area 0{u] r). When the reference point 



is located at the center of the L-gon, the overlap area 0{u]r) can be easily evaluated as shown in [[T|. 
For < r < i?i (where Ri is the L-gon inradius defined above ([T])), the disk V{u] r) is entirely inside the 
L-gon A. Thus the overlap area is the same as the area of the disk, i.e, 0{u] r) = vrr^. For Ri < r < R, 
(where R is the L-gon circumradius) there are L circular segment shaped portions of the disk V{u]r) 
which are outside the L-gon. Since these circular segments are symmetric and there is no overlap between 
them, the overlap area is given by 0{u]r) = vrr^ — L(area of one circular segment). Note that in this 
case it is simple to find the area of one circular segment using geometry. 

When the reference point is not located at the center, then the circular segments are no longer symmetric. 
This is illustrated in Fig.[2]for the case of a square (L = 4) for simplicity. We can see that for the given 
radius of the disk V{u] r) centered at the reference point, non-symmetric circular segment areas are formed 
outside the sides Si, S2 and 6^4. In addition, there is overlap between the circular segment areas formed 
outside the sides Si and 5^2. This complicates the task of finding the overlap area 0{u]r). Also, for any 
arbitrary location of the reference point, the different ranges for r are no longer exclusively determined 
by the inradius Ri and circumradius i?, but by the number of unique distances from the reference point 
to the vertices and the sides. This adds further complexity to the task of finding the overlap area 0{u] r). 

In the next section, we propose our method to systematically account for the effect of sides and vertices 
on the overlap area 0{u]r). Then in Section |v| we propose an algorithm to automatically formulate the 
different ranges for r and find the overlap area 0{u] r) in order to evaluate ( [TT] ) and ( [T2| ). 



IV. Characterization of the Effect of Sides and Vertices 

In this section, we characterize the effect of polygon sides and vertices on the overlap area 0{u]r). 
Because of the symmetry of the polygon, we only need to consider the following three cases which are 
illustrated in Fig.[3| Other cases can be handled as an appropriate combination of these cases: 

1) The overlap area is limited by one side of the polygon only. This is illustrated in Fig.[3]^a) for side 

2) The overlap area is limited by two sides of the polygon and there is no overlap between them. This 
is illustrated in Fig.|3]^b) for sides Si and Sl- 

3) The overlap area is limited by two sides of the polygon and there is overlap between them. This is 
illustrated in Fig.[3]^c) for side Si and Sl and vertex Vi. 

These cases are discussed in detail in the following subsections. 

A. Case 1: Overlap area limited by one side only 

Let Bi{u] r) denote the circular segment formed outside the side Si, as shown in Fig.[3]^a). It is obvious 
that the overlap area in this case is given by 0{u]r) = vrr^ — Bi{u]r), where Bi{u]r) = \Bi{u]r)\, 
Using polar coordinates, we can find the area Bi{u]r) by integrating the angle 6i (shown in Fig.jsj^a)) 
over Bi{u] r) as 



Bi{u]r) = \\Bi{u]r) 



= 2 / arccos [ ^^^^^^^^ \ ^j^' 

^2a,,cos (^^) - {d{u; S^)f arccos - 

5i) (^Y^r2 - {p{u; S^))' - ^ S^)y - (p(ti; , r>diu;S,); (13) 

0, otherwise; 
where p{u] Si) and d{u] Si) are given in ^ and ([8]) respectively. 

Generalizing, let B£{u] r) denote the circular segment which is formed outside the side S^. Using ([13]) 
and the rotation operator in ([3]), we can express B^^u; r) as 

5i(9^-(^-i)i6;r), r > d{u]Si)] 
B,{u;r) = { (14) 

0, otherwise. 



= < 



B. Case 2: Overlap area limited by two sides with no overlap 

This case is illustrated in Fig.jsj^b), where two circular segments, Bi{u]r) and BL{u]r), are formed 
outside the sides Si and Sl, respectively. Since there is no overlap between the circular segments, the 
overlap area is given by 0{u]r) = vrr^ — {Bi{u]r) + Bl^u^t)), where Bi{u]r) is given by ( [T3] ) and 
BL{u]r) can be found using ([14]). 



C. Case 3: Overlap area limited by two sides with overlap 

This case is illustrated in in Fig.[3]^c), where there is overlap between the two circular segments, Bi{u] r) 
and BL{u]r), formed outside sides Si and Sl due to the inclusion of the vertex Vi. Let Ci{u]r) = 
Bi{u] r) n Bl{u] r) denote this overlap region between the two circular segments. Thus, in this case, the 
overlap area is given by 0{u] r) = vrr^ — {Bi{u] r) + Bl{u] r) — Ci{u] r)). 

Using polar coordinates, we can find the area Ci{u] r) by integrating the angle 62 (shown in Fig.[3]^c)) 
over Ci{u] r) as 



Ci(w;r) = ||Ci(w;r) 



d{u-Vi) 



L 



+ arccos 



p{u]Si) 



+ arccos 



p{u]Sl) 



TT 1 dr' 



— I arccos ( I ' ) + arccos ' ' 



(arccos (gfj 



2 



arccos 



+ 



[^j{d{u', Vi))" - {p{u; Si)Y - - {p{u; Si)Y^ 
^ - {d{u;Vi)Y) ^ 



(15) 



^ y^y^, , , r> d{u]Vi)] 

0, otherwise; 
where d{u] Vi), p{u] Si) and p{u] Sl) are given in ([5]), ^ and ^ respectively. 

Generalizing, let Ci{u] r) denote the overlap region between two circular segments adjoining the vertex 
Vi of the polygon. Using ([15]) and the rotation operator in (|3]), we can express Ci{u] r) as 

Ci{pK-^^-^^u',r), r > d{u]Vi)] 

(16) 

0, otherwise. 



Cpiu\r) 



The expressions for the derivatives of Bi{u] r) and Ci{u] r), which are needed in the evaluation of ^^^^'^^ 
in ( [T2| ), are provided in Appendix|A[ 



V. Distance Distributions in Polygons 



In this section, we present our algorithm to use the distance vector in ([TO]) and the effect of different 
sides and vertices ((p3])— ([T6|)) in order to find the overlap area 0{u\r) given any arbitrary reference point 
u = [x, nY located inside the L-gon A, 

The overall effect of the sides and vertices of the L-gon depends on the range of the distance r and 
on the distance between the reference point and all the sides and vertices. Since an L-gon has L sides 
and L vertices, there can be 2L + 1 different ranges for the distance r. We sort the distance vector d in 



([To]) in ascending order and define d to be the sorted distance vector. Then, the first range of the distance 
is < r < (ii, where di denotes the first entry of the sorted distance vector d. The next 2L — 1 ranges 
are dj < r < rfj+i, j = 2, 3, . . . , 2L and the last range is d2L < Thus, in general, we can write an 
expression for 0{u]r) as 

Oi{u] r) = 7rr^, < r < rfi; 

02{u]r), di <r < 

0{u]r) = { \ : (17) 

02l{u] r), d2L-i <r < ^l; 

02L+i{u]r) = A, r> rfsL; 
where the subscript j in Oj{u] r) denotes the overlap area for the j-th range. 

For the first range Q < r < di, Si{u]r) = vrr^ since the disk V{u]r) will be completely inside the 
^-gon. For the last range d2L < 02L+i{u]r) = A since the disk V{u]r) will completely overlap with 
the L-gon A. For the intermediate ranges, the overlap area may be limited by any number of sides, with or 
without overlap between any two adjacent sides. Also depending on the location of the arbitrary reference 
point, some of the distance to the sides and vertices may be the same. In order to automate the process 
of finding the unique set of ranges for the distance r and to calculate the corresponding overlap area for 
each unique range, we propose the Algorithm [Tj 



In the proposed algorithm, we sort the distance vector d in ( [TO] ) in ascending order to obtain the sorted 



distance vector d and the index vector k = [fci, /c2, . . . , /c2l]. The index vector k is then used to find the 
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Algorithm 1 Algorithm to find the overlap area 



Step 1: Sort d in ( [T0| ) in ascending order to obtain d 

Step 2: Determine the index sorting that transforms d into d and obtain the index vector k 
Step 3: Find the appropriate circular segment areas and the overlap area 
for each j in j = 1, 2, 3, . . . , 2L + 1 do 
if dj-i — dj ^ 0, (rfo = 0) then 
0^{u\ t) = 7rr^ 

for each i in i = 1, 2, 3, . . . , j — 1 do 
\iki<L then 

Oj{u',r) = Oj{u]r) - B^Xu^r) 
else 

Oj{u]r) = Oj{u]r) + Ck,-L{u]r) 
end if 
end for 
end if 
end for 



effect of sides and vertices on the overlap area for each value of range. For each unique range dj^i — dj 



which is determined in the outer for loop, we evaluate the overlap area Oj{u] r) in ( ^7| ) in the inner for 
loop by taking into account the cumulative effect of all of the sides and vertices. The terms Bl{u] r) and 
Cl{u] r) are appropriately selected using the index vector k. Once 0{u] r) in ( [T7| ) is computed for each 
unique range, it can then used to evaluate ( pTj ) and ([12]). We have implemented the proposed algorithm 
in MATLAB. The code is provided as supplementary material with this paper. 



A. Algorithm illustration - arbitrary reference point in a square 

In order to illustrate the proposed framework and algorithm, we consider the problem to find CDF in 
( [TT] ) and the PDF in ([12]) for a point Ui = [R/2, -R/2f located in the middle of the side S4 of the 



square (L = 4). The distance vector d in ( fTO) ) and the sorted distance vector d are given by 

VTOR R_ 



R R ^ R 



JL JL JL JL /^p yio-^ vioi^ 
' 71' 71' 71' 71' 



Thus the index vector k for this case is given by k = [4, 1, 2, 5, 8, 3, 6, 7], which is employed to 



determine 0{u]r) in ( [17] ) using the proposed algorithm. Thus the CDF is 



7rr 



(S1 + E2 + 54-C1-C4), 



< r < i?/y2; 



(18) 



7rr2 - (5i + ^2 + ^3 + 54 - Ci - C4), ^i? < r < ^i?/2; 



r > V10i?/2. 



By plugging the CDF F{u;r) in ([TS]), and its derivative which is obtained using the expressions in 
Appendix|A] in ([12]), we obtain the PDF of the distance to the n-th neighbor, which is shown in Fig.|4] 
for i? = 1, = 5 and n = 1, 2,..., 5. We can see that the simulation results, which are averaged over 
10,000 simulation runs, match perfectly with the analytical results, verifying the proposed framework and 
algorithm. 



B. Special case - reference point at the center of the polygon 

Consider the special case that the reference point is located at the center O of the polygon, i.e., 
at the origin U2 = [0 0]^. All of the sides and also all of vertices are equidistant from the center 
of the polygon and the rotation operator does not affect the point located at the origin. This implies 
that d{ui]S^) = d{ui;Si) = Rtan{e/2), d{ui;Ve) = d{ui;Vi) = R, Be = Bi and Cg = Ci for 
^ = 2, 3, . . . , L. By employing these relations and the proposed algorithm, we obtain the CDF expression 
for the following three possible ranges 



vrr 



F{u2;r) = - <j 



< r < i?tan(|) ; 
7rr2 - LBi, i?tan (f ) < r < R; 
A, r> R; 



(19) 



where Bi = arccos tan (|) ) — i?tan | 
in ( fT2l ) reproduces the result in |[1|. 



i?2tan2 (I) and 6 is defined in Substituting ([19]) 
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C. Special case - reference point at the vertex of the polygon 

Consider the special case when the arbitrary point is located at one of the vertex of the polygon. Let 
the point be located at Vi, that is = [R 0]^. From the vertex Vi, the sides S^, Sl+i-i and the 
vertices V^, are all equidistant for ^ = 2, 3, . . . [L/2J. Hence, there are | + 1 possible ranges of 

the distance. By using the proposed algorithm to determine the overlap area 0{u] r), we obtain the CDF 
expression for these ranges as 

i=l 



d{u,, V,)<r< d{u,, ^ = 1, 2, . . . Lf J , Cl+, = 0; (20) 



which is then used to the evaluate the PDF fn{u] r) of the distance to the n-th neighbor by employing ( [T2| ). 
Fig. [5] shows the plot of the PDF of the distance to the farthest neighbour (n = 10) from the vertex Vi, 
with = 10 nodes distributed inside a ^-gon with area A = 100, for ^ = 3,4, 6, 100. For the simulation 
results, we have used the method in [15J to uniformly distribute the points inside the ^-gon and averaged 
the results over 10,000 simulation runs. Again it can be seen that the simulation results match perfectly 
with the analytical results. 

VI. Conclusions 

In this paper, we have proposed a general framework for analytically obtaining the closed-form proba- 
bility density function of the distance between any arbitrary reference point and its n-th neighbour node, 
when nodes are uniformly and independently distributed inside a regular L-sided polygon. We have 
provided examples to demonstrate the generality of our proposed framework. It is noted that the general 
distance distributions obtained in this paper can be applied to accurately model and analyse the wireless 
network characteristics such as connectivity, interference and outage probability. 
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Appendix A 
Derivatives 

By employing the Leibniz integral rule, we can express the derivatives of Bi{u] r) and Ci{u] r), which 
are required in the evaluation of ^^'^'^'^ in ( [T2| ), as 

^^'^^-''■^ -2r.rccoJ'-^) (2.) 



dr \ r 

dCe{u-r) ( [ p{u-S,) \ , (p{u;^^\ 2n 

rl arccos ^ + arccos ^ — (22) 



dr V V^/ \ T J L 



References 

[1] S. Srinivasa and M. Haenggi, "Distance distributions in finite uniformly random networks: Theory and applications," IEEE Trans. Veh. 

Technol, vol. 59, no. 2, pp. 940-949, Feb. 2010. 
[2] D. Moltchanov, "Distribution of distances in random networks," Ad Hoc Networks, Mar. 2012. [Online]. Available: 

|http://dx.doi.org/10.1016/j.adhoc.2012.02.005| 
[3] P. Fan, G. Li, K. Cai, and K. B. Letaief, "On the geometrical characteristic of wireless ad-hoc networks and its application in network 

performance analysis," IEEE Trans. Wireless Commun., vol. 6, no. 4, pp. 1256-1265, Apr. 2007. 
[4] M. Haenggi and R. K. Ganti, Interference in Large Wireless Networks. Now Publishers Inc., 2008. 
[5] J. G. Andrews and S. Webber, Transmission capacity of wireless networks. Now Publishers Inc., 2012. 

[6] F. Fabbri and R. Verdone, "A statistical model for the connectivity of nodes in a multi-sink wireless sensor network over a bounded 

region," in Proc. 14th European Wireless Conference, Jun. 2008. 
[7] M. Haenggi, "On distances in uniformly random networks," IEEE Trans. Inf. Theory, vol. 51, no. 10, pp. 3584-3586, Oct. 2005. 
[8] D. Torrieri and M. C. Valenti, "The outage probability of a finite ad hoc network in Nakagami fading," IEEE Trans. Commun., 2012 

(to appear). [Online]. Available: [http://arxiv.org/abs/1207.271 lvl| 
[9] G. Mao and B. D. O. Anderson, "Towards a better understanding of large-scale network models," IEEE/ACM Trans. Netw., vol. 20, 

no. 2, pp. 408-421, Apr. 2012. 

[10] J. Guo and S. Durrani, "Probability distribution of distance between wireless nodes in a triangular region," Electronics Letters, Apr. 
2012 (under review). 

[11] A. M. Mathai, An introduction to geometrical probability: distributional aspects with applications. CRC Press, 1999. 
[12] A. Kostin, "Probability distribution of distance between pairs of nearest stations in wireless network," Electronics Letters, vol. 46, 
no. 18, pp. 1299-1300, Sep. 2010. 

[13] Y. Zhuang and J. Pan, "Random distances associated with hexagons," ArXiv Technical Report, 2011. [Online]. Available: 
|http://arxiv.org/abs/l 106.2200| 

[14] U. Basel, "Random cords and point distances in regular polygons," ArXiv Technical Report, 2012. [Online]. Available: 
|http://arxiv.org/abs/1204.2707| 

[15] M. Hlynka and D. Loach, "Generating uniform random points in a regular n sided polygon," Windsor Math and Statistics Report, 
2005. [Online]. Available: |http://www2.uwindsor.ca/^hlynka/loachpolygoiLpdf, 



12 



List of Figures 

1 Illustration of sides, vertices and angles for the L-sided polygon inscribed in a circle of radius 
R. 6 denoting the interior angle of the polygon and i9 denoting the central angle between 

two adjacent vertices are defined in ([2]) [13] 

2 Illustration of the non- symmetric circular segment areas formed outside the sides Si, S2 and 
5^4 for a square (L = 4-gon). There is also an overlap between the circular segment areas 
formed outside the sides Si and S2 [H 

3 Illustration of the effect of sides and vertices of the polygon on the overlap are: (a) overlap 
area is limited by side Si, (b) overlap area is limited by sides Si and Sl only and (c) overlap 

area is limited by sides Si and Sl with inclusion of the vertex Vi [15] 

4 The PDF fn{u] r) of the distance to n-th nearest neighbor for a point located midway between 
the Vi and V4 on the side 6^4 of the square (L = 4-gon). The circumradius i? = 1 and the 
number of nodes is = 5 [H] 

5 The PDF /io(i^;r) of the distance to farthest neighbor for a point located at the vertex 
of the polygon with area A = 100 and = 10 nodes, for ^ = 3,4, 6, 100. Solid lines 
show analytical result using our proposed algorithm and markers provide verification using 
simulation [17] 
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Fig. 1. Illustration of sides, vertices and angles for the L-sided polygon inscribed in a circle of radius R. 6 denoting the interior angle of 
the polygon and i9 denoting the central angle between two adjacent vertices are defined in (|2|. 
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Fig. 2. Illustration of the non- symmetric circular segment areas formed outside the sides Si, S2 and *S'4 for a square (L = 4-gon). There 
is also an overlap between the circular segment areas formed outside the sides 5*1 and 6*2 . 
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(c) 

Fig. 3. Illustration of the effect of sides and vertices of the polygon on the overlap are: (a) overlap area is limited by side Si, (b) overlap 
area is limited by sides Si and Sl only and (c) overlap area is limited by sides 5*1 and Sl with inclusion of the vertex Vi. 
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3.5 




distance, r 

Fig. 4. The PDF fn{u; r) of the distance to n-th nearest neighbor for a point located midway between the Vi and V4 on the side & of 
the square (L = 4-gon). The circumradius R — 1 and the number of nodes is N = 5. 
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Fig. 5. The PDF fio{u;r) of the distance to farthest neighbor for a point located at the vertex of the polygon with area A — 100 and 
= 10 nodes, for ^ = 3, 4, 6, 100. Solid lines show analytical result using our proposed algorithm and markers provide verification using 
simulation. 



